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Preface 



This book is not meant to be another compendium of select inequalities, nor does it claim 
to contain the latest or the slickest ways of proving them. This project is rather an attempt at 
describing how most functional inequalities are not merely the byproduct of ingenious guess 
work by a few wizards among us, but are often manifestations of certain natural mathematical 
structures and physical phenomena. Our main goal here is to show how this point of view 
leads to "systematic" approaches for not just proving the most basic functional inequalities, 
but also for understanding and improving them, and for devising new ones - sometimes at 
will, and often on demand. 

Our aim is therefore to describe how a few general principles are behind the validity of 
large classes of functional inequalities, old and new. As such, Hardy and Hardy-Rellich type 
inequalities involving radially symmetric weights are variational manifestations of Sturm's the- 
ory on the oscillatory behavior of certain ordinary differential equations. Similarly, allowable 
non-radial weights in Hardy-type inequalities for more general uniformly elliptic operators 
are closely related to the resolution of certain linear PDEs in divergence form with either a 
prescribed boundary condition or with prescribed singularities in the interior of the domain. 

On the other hand, most geometric inequalities including those of Sobolev and Log- 
Sobolev type, are simply expressions of the convexity of certain free energy functionals along 
the geodesies of the space of probability measures equipped with the optimal mass transport 
(Wasserstein) metric. Hardy-Sobolev and Hardy-Rellich-Sobolev type inequalities are then 
obtained by interpolating the above inequalities via the classical ones of Holder. 

Besides leading to new and improved inequalities, these general principles offer new ways 
for estimating their best constants, and for deciding whether they are attained or not in 
the appropriate function space. In Hardy-type inequalities, the best constants are related to 
the largest parameters for which certain linear ODEs have non-oscillatory solutions. Duality 
methods, which naturally appear in the new "geodesic convexity" approach to geometric 
inequalities, allow for the evaluation of the best constants from first order equations via the 
limiting case of Legendre-Fenchel duality, as opposed to the standard method of solving second 
order Euler-Lagrange equations. 

Whether a "best constant" on specific domains is attained or not, is often dependent on 
how it compares to related best constants on limiting domains, such as the whole space or 
on half-space. These results are based on delicate blow-up analysis, and are reminiscent of 
the prescribed curvature problems initiated by Yamabe and Nirenberg. The exceptional case 
of the Sobolev inequalities in two dimensions initiated by Trudinger and Moser can also be 
proved via mass transport methods, and some of their recent improvements by Onofri, Aubin 
and others are both interesting and still challenging. They will be described in the last part 
of the monograph. 

The part dealing with Hardy and Hardy-type inequalities represents a compendium of 
work mostly done by -and sometimes with- my (now former) students Amir Moradifam and 
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Craig Cowan, while the "mass transport" approach to geometric inequalities follows closely 
my work with my former student X. Kang and postdoctoral fellow Martial Agueh. This is 
largely based on the pioneering work of Cedric Villani, Felix Otto, Robert McCann, Wilfrid 
Gangbo, Dario Cordero-Erausquin, B. Nazareth, C. Houdre and many others. The chapters 
dealing with Hardy-Sobolev type inequalities follow work done with my students Chaogui 
Yuan, and Xiaosong Kang, a well as my collaborator Frederic Robert. Finally, much of the 
progress on the -still unresolved- best constant in Moser-Onofri-Aubin inequalities on the 
2-dimensional sphere was done with my friends and collaborators, Joel Feldman, Richard 
Froese, Changfcng Gui, and Chang-Shou Lin. I owe all these people a great deal of gratitude. 

I thank my wife Louise and my children Mireille, Michelle and Joseph for their patience 
and their support over the years. 



Nassif Ghoussoub 



Introduction 



This book is an attempt to describe how a few general principles are behind the validity 
of large classes of functional inequalities, old and new. It consists of six parts, which -though 
interrelated- are meant to reflect either the mathematical structure or the physical phenomena 
behind certain collections of inequalities. 

In Part I, we deal with Hardy-type inequalities involving radially symmetric weights and 
their improvements. The classical Hardy inequality asserts that for a domain f2 in K™, n > 3, 
with £ fi, the following holds: 

(0.1) f n \Vu\Hx>{^Yj n ^dx tovu£Hl{n). 

The story here is the newly discovered link between various improvements of this inequality 
confined to bounded domains and Sturm's theory regarding the oscillatory behavior of certain 
linear ordinary equations, which we review in Chapter 1. 

In Chapter 2, we first identify suitable conditions on a non-negative C 1 -function P defined 
on an interval (0, R) that will allow for the following improved Hardy inequality to hold on 
every domain ft contained in a ball of radius R: 

(0.2) /„ {^f f a ^dx > J n P(\x\)u 2 dx for u £ H^Q). 

It turned out that a necessary and sufficient condition for P to be a Hardy Improving Potential 
(abbreviated as Hl-potentiat) on a ball Br, is for the following ordinary differential equation 
associated to P 

(0.3) y" + V + P(r)y = 0, 

r 

to have a positive solution on the interval (0,R). Elementary examples of Hl-potentials are 
P = on any interval (0,R), P = 1 on (0,zo), where zq = 2.4048... is the first zero of the 
Bcssel function J , and more generally P(r) = r~ a with < a < 2 on (0, z a ), where z a is the 
first root of the largest solution of the equation y" + -y' + r~ a y = 0. Other examples are 

P P (r) - ir ^ ogtyl on (0, f ), but also P fc » = 4, f (UU log® f)~ 2 on (0, e( /_ Umes) ). 

Besides leading to a large supply of explicit Hardy improving potentials, this connection 
to the oscillatory theory of ODEs, gives a new way of characterizing and computing best 
possible constants such as 



(0.4) P(P,R):= mf 



f | Vu | 2 dx-^—^- ( \x\- 2 \u\ 2 dx 



ueH^n) J Q P(\x\)u 2 dx 

M#0 
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On the other hand, the value of the following best constant 

/ |Vu| 2 dx-X f P(\x\)u 2 dx 
Jn Jn 



(0.5) nx(P,n):= inf 



Jn 



and whether it is attained, depend closely on the position of the singularity point vis-a-vis 
0. It is actually equal to ^ n ~^ , and is never attained in Hq(Q), whenever fi contains in 
its interior, but the story is quite different for domains f2 having on their boundary. In this 
case, fj,\(P,£i) is attained in Hq(Q.) whenever ^\(P,fl) < which may hold or not. For 

example, fi\(P, 0) is equal to ^ for domains that lie on one side of a half-space. 

In Chapter 3, we consider conditions on a couple of positive functions V and W on (0, oo), 
which ensure that on some ball Br of radius R in K", n > 1, the following inequality holds: 

(0.6) J B V(\x\)\\7u\ 2 dx > J B W(\x\)u 2 dx for u e C^{B R ). 

A necessary and sufficient condition is that the couple (V, W) forms a n- dimensional Bessel 
pair on the interval (0, R), meaning that the equation 

/ n ///x ,n—l V r (r) . ., W(r) , , 

(0.7) y"(r) + (— - + y^)y'(r) + ^y(r) = 0, 

has a positive solution on (0, R). This characterization allows us to improve, extend, and unify 
many results about weighted Hardy-type inequalities and their corresponding best constants. 
The connection with Chapter 2 stems from the fact that P is a Hi-potential if and only if the 
couple (1, + P) is a Bessel pair. More generally, the pair 

(0.8) ( r-\ C~^~ 2 ) 2 r- X - 2 + r- x P(r) 



is also a n-dimensional Bessel pair on (0, R) provided < A < n — 2. Again, the link to Sturm 
theory provides many more examples of Bessel pairs. 

Hardy's inequality and its various improvements have been used in many contexts such as 
in the study of the stability of solutions of semi-linear elliptic and parabolic equations, of the 
asymptotic behavior of the heat equation with singular potentials, as well as in the stability 
of eigenvalues for Schrodinger operators. In Chapter 4, we focus on applications to second 
order nonlinear elliptic eigenvalue problems such as 



(0.9) 



— Au — Xf(u) in Q 
u = on dfl, 



where A > is a parameter, f2 is a bounded domain in WL N , N > 2, and / is a superlinear 
convex nonlinearity. The bifurcation diagram generally depends on the regularity of the 
extremal solution, i.e., the one corresponding to the largest parameter for which the equation 
is solvable. Whether, for a given nonlinearity /, this solution is regular or singular depends 
on the dimension, and Hardy-type inequalities are crucial for the identification of the critical 
dimension. 

Part II deals with the Hardy-Rellich inequalities, which are the fourth order counterpart 
of Hardy's. In Chapter 5, we show that the same condition on the couple (V, W) (i.e, being a 
n-dimensional Bessel pair) is also key to improved Hardy-Rellich inequalities of the following 
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type: For any radial function u G C§°(Br) where Br is a ball of radius R in R n , n > 1, we 
have 

(0.10) J B V(\x\)\Au\*dx > f B W(\x\)\Vu\ 2 dx + (n - 1) J B (^ - ^)|V U | 2 dz. 
Moreover, if 

(0.11) W(r)-^ + ^-V rr (r)>0 on[0,R), 

then the above inequality holds true for all u € C§ c (Br) and not just the radial ones. By 
combining this with the inequalities involving the Dirichlet integrals of Chapter 3, one ob- 
tains various improvements of the Hardy- Rellich inequality for ffg(O). In particular, for 
any bounded domain Cl containing with 57 C -Br, we have the following inequality for all 
u€H 2 (Q), 



(0 , 2) 

where n > 4, A < n — 2, and where P is a Hi-potential on (0, R) such that ^£jy = 7 + /M, 
/(r) > and lim rf(r) = 0. 

In Chapter 6, we explore Hardy-type inequalities for iJ 1 (ri)-functions, i.e., for functions 
which do not necessarily have compact support in f2. In this case, a penalizing term appears 
in order to account for the boundary contribution. If a pair of positive radial functions (V, W) 
is a n-dimensional Bessel pair on an interval (0,i?), and if Br is a ball of radius R in R™, 
n > 1, then there exists 9 > such that the following inequality holds: 

(0.13) J Br V(x)\Vu\ 2 dx > J Br W{x)u 2 dx - 6 J 9Br u 2 ds for u 6 H l {B R ), 
and for all radial functions u £ H 2 (Br), 

(0.14) / V(\x\)\Au\ 2 dx > f W(\x\)\Vu\ 2 dx + (n-l) f _ M^l)| Vu | 2 dx 

Jb r Jb r iB R pi Fl 

+ [(n-l)-0)V(J?)] / |Vw| 2 da;. 



SB/ 



The latter inequality holds for all functions in H 2 (B) provided (|0.11[) holds. The combination 
of the two inequalities lead to various weighted Hardy- Rellich inequalities on H 2 n Hq . 

In Chapter 7, we investigate some applications of the improved Hardy- Rellich inequalities 
to fourth order nonlinear elliptic eigenvalue problems of the form 



(0.15) 



A 2 u = Xf(u) in fl 
u = Au = on dfl, 



as well as their counterpart with Dirichlet boundary conditions. In particular, they are again 
crucial for the identification of "critical dimensions" for such equations involving either an 
exponential or a singular supercritical nonlinearity. 

Part III addresses Hardy-type inequalities for more general uniformly elliptic operators. 
The issue of allowable non-radial weights (to replace j^p) is then closely related to the res- 
olution of certain linear PDEs in divergence form with either prescribed conditions on the 
boundary or with prescribed singularity in the interior. We also include L p -analogs of various 
Hardy- type inequalities. 
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In Chapter 8, the following general Hardy inequality is associated to any given symmetric, 
uniformly positive definite nxn matrix A{x) defined in SI with the notation := (A(x)£, £) 
for £ G M™. 

(0.16) / \Vu\\dx >- 1 l A u 2 dx, u e i^(fi) 

Jn 4 J n E 

The basic assumption here is that E is a positive solution to — div(AVE) dx — fi on Jl, where 
/i is any nonnegative nonzero finite measure on fl. The above inequality is then optimal in 
either one of the following two cases: 

• E is an interior weight, that is E — +oo on the support of fj,, or 

• E is a boundary weight, meaning that E — on dQ. 

While the case of an interior weight extends the classical Hardy inequality, the case of a 
boundary weight extends the following so-called Hardy 's boundary inequality, which holds for 
any bounded convex domain Jl C K™ with smooth boundary: 

(0.17) J n \Vu\ 2 dx > \ J Q dist £ dn)2 dx for u e ff i(fi). 

Moreover the constant j is optimal and not attained. One also obtains other Hardy inequal- 
ities involving more general distance functions. For example, if Jl is a domain in R" and M 
a piecewise smooth surface of co-dimension k (k — l,...,n). Setting d(x) := dist(x, M) and 
suppose k ^ 2 and ~Ad 2 ~ k > in Cl\M, then 

(0.18) /„ \Vu\ 2 dx > J n ^dx for u G ^(flXM). 

The inequality is not attained in either case, and one can therefore get the following improve- 
ment for the case of a boundary weight: 

(0.19) f \Vu\\dx > - [ l -^^u 2 dx + l -f ^d/i, u G Hl(n) 

which is optimal and still not attained. Optimal weighted versions of these inequalities are 
also established, as well as their L p -counterparts when p ^ 2. Many of the Hardy inequalities 
obtained in the previous chapters can be recovered via the above approach, by using suitable 
choices for E and A(x). 

In Chapter 9, we investigate the possibility of improving (|0.16[) in the spirit of Chapters 
3 and 4, namely whether one can find conditions on non-negative potentials V so that the 
following improved inequality holds: 

(0.20) / n \Vu\ 2 A dx ~y Q ^^u 2 dx > J Q V{x)u 2 dx for u G 

Necessary and sufficient conditions on V are given for (|0.20l) to hold, in terms of the solvability 
of a corresponding linear PDE. Analogous results involving improvements are obtained for 
the weighted versions. Optimal inequalities are also obtained for i? 1 (S7). 

We conclude Part III by considering in Chapter 10, applications of the Hardy inequality 
for general uniformly elliptic operators to study the regularity of stable solutions of certain 
nonlinear eigenvalue problems involving advection such as 

—Am + c(x) ■ Vu = (i- u )2 m 
u = on dSl, 



(0.21) 



where c(x) is a smooth bounded vector field on 51. 
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In Part IV, we describe how the Monge-Kantorovich theory of mass transport provides 
a framework that encompasses most geometric inequalities. Of importance is the concept of 
relative energy of p with respect to p\ defined as: 



(0-22) Rp w (p \p 1 ):=Up w (p )-Rp w (p 1 ), 

where p and p\ are two probability densities, and where the Free Energy Functional Hy W 
is defined on the set V a {£l) of probability densities on a domain f2 as: 



(0.23) T# w {p) := [ 

Jn 



F{p)+ P V+ l -(W*p)p 



Ax. 



~Ry W being the sum of the internal energy H F (p) := J n F(p)dx, the potential energy Hy (p) := 
J Q pVdx and the interaction energy K w (p) := | f Q p(W * p) dx. Here F is a differentiable 
function on (0, oo), while the confinement (resp., interactive) potential V (resp., W) are C 2 - 
functions on R n satisfying D 2 V > pi (resp., D 2 W > vT) for some p, v e M. 

In Chapter 11, we describe Brenier's solution of the Monge problem with quadratic cost, 
which yields that the Wasserstein distance W(po,pi) between two probability densities po on 
X and p\ on Y, i.e., 

(0.24) W(p , Pl ) 2 =m{^J^\x-s{x)\ 2 dx; sGS(po,Pi)} 

is achieved by the gradient V<p of a convex function tp. Here S(po, p\) is the class of all Borel 
measurable maps s : X — > Y that "push" p into pi, i.e., those which satisfy the change of 
variables formula, 

(0.25) / h(y)p 1 (y)dy= / h(s(x))p (x)dx for every h G C(Y). 

This fundamental result allows one to show that for certain natural candidates F, V and 
W, the corresponding free energy functionals Hy w are convex on the geodesies of optimal 
mass transport joining two probability densities in V a {Q)- This convexity property translates 
into a very general inequality relating the relative total energy between the initial and final 
configurations po and p\ 1 to their entropy production I c *(p\p v ), their Wasserstein distance 
W 2 (po, pi), as well as the Euclidean distance between their barycenters |b(po) — b(yOi)|, 



(0.26) H^(a,| Pi ) + ^W 2 (po, Pl ) £|b(po) -b( Pl )| 2 < KlvvT W (Po) + I c *(p\p v ). 



Here Pf{x) := xF'(x) — F(x) is the pressure function associated to F, while c is a Young 
function (such as c(x) — ||a;| p ), c* is its Legendre transform, while I c *(p\p v ) is the relative 
entropy production-type function of p measured against c* defined as 

(0-27) I c *{p\p v ):= [ pc*(-\7(F'(p) + V + W*p))dx. 

Once this general comparison principle is established, various - new and old - inequalities 
follow by simply considering different examples of internal, potential and interactive energies, 
such as F(p) = p\np or F(p) = p 1 , and V and W are convex functions (e.g., V(x) = jM 2 ), 
while W is required to be even. 

The framework is remarkably encompassing even when V — W = 0, as it is shown in 
Chapter 12 that the following inequality, which relates the internal energy of a probability 
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density p on R™ to the corresponding entropy production contains almost all known Euclidean 
Sobolev and log-Sobolev inequalities: 

(0.28) f [F(p) + nP F (p)]dx < [ pc*(-V(F'op)) dx + K c . 

Jn Jn 

The latter constant K c can always be computed from F and the Young function c. 

The approach allows for a direct and unified way for computing best constants and ex- 
tremals. It also leads to remarkable duality formulae, such as the following associated to the 
standard Sobolev inequality for n > 3 and where 2* := ^35: 

(0.29) snp{^—p-[ p{x)^dx- [ \x\ 2 p(x)dx; [ p{x) dx = 1} 

= inf{ f \\7f\ 2 dx; f 6 C* °°(R"); f \ffdx = l}. 

L JR™ JR" ' 

This type of duality also yields a remarkable correspondence between ground state solutions 
of certain quasilinear (or semi-linear) equations, such as "Yamabe's" , 

-A/ = |/| 2 *- 2 / on R", 

and stationary solutions of the (non-linear) Fokker-Planck equations — Au 1_ ™ + div(x.it), 
which -after appropriate scaling- reduces to the fast diffusion equation 

= Au 1 -^ on R+ x R n . 

dt 

Chapter 13 deals with applications to Gaussian geometric inequalities. We first establish 
the so-called HWBI inequality, which follows immediately from a direct application of (|0.26|) 
with paramatrized quadratic Young functions c a (x) — ^\%\ 2 for a > 0, coupled with a 
simple scaling argument: 

(0.30) R$ W (p \pi) < W 2 (p , Pl Wl 2 (p \p v ) a±^wi(po, Pl ) + ||b(p ) -b( Pl )| 2 . 

This gives a unified approach for -extensions of- various powerful inequalities by Gross, 
Bakry-Emery, Talagrand, Otto-Villani, Cordero-Erausquin, and others. As expected, such 
inequalities also lead to exponential rates of convergence to equilibria for solutions of Fokker- 
Planck and McKean-Vlasov type equations. 

Part V deals with Caffarelli-Kohn-Nirenberg and Hardy-Rellich-Sobolev type inequalities. 
All these can be obtained by simply interpolating -via Holder's inequalities- many of the 
previously obtained inequalities. This is done in Chapter 14, where it is also shown that the 
best constant in the Hardy-Sobolev inequality, i.e., 

(0.31) /i s (f2) := inf / f |Vu| 2 g?:e; u G Hq(Q) and / ^ / } dx = l| , 

Un Jn FT J 

where < s < 2 and 2*(s) — 2 ( n-s ) j s never attained when is in the interior of the domain 
f2, unless the latter is the whole space R™, in which case explicit extremals are given. This is 
not the case when ft is half-space R", where only the symmetry of the extremals is shown. 
Much less is known about the extremals in the Hardy-Rellich-Sobolev inequality (i.e., when 
s > 0) even when £1 = R n . 

The problem whether p s (£l) is attained becomes more interesting when is on the bound- 
ary d£l of the domain il. The attainability is then closely related to the geometry of dCl, as 
we show in chapter 15, that the negativity of the mean curvature of dfl at is sufficient to 
ensure the attainability of p s (Q). 
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In Chapter 16, we consider log-Sobolev inequalities on the line, such as those involving 
the functional 

(0.32) I a (g) = | J^(l-x 2 )\g'(x)\ 2 dx + j g{x) dx-\n\ j e 2 ^ dx 

on the space 1, 1) of L 2 -functions on (—1, 1) such that — x 2 )\g' (x)\ 2 dx) x / 2 < oo. 

We then show that if J a is restricted to the manifold 

G = {<? € 1); J 1 e 29 ^xdx = o| , 

then the following hold: 

(0.33) inf I a (g) = if a > \, and inf I a (g) = — oo if a < ^ 

We also give a recent result of Ghigi, which says that the functional 
(0.34) $( u )=y u(x) dx - \og(J^ J + e- 2u * {x Ux^J 

is convex on the cone W of all bounded convex functions u on (— 1, 1), where here u* denotes 
the Legendre transform of u, and that 

4 

inf : $(«) = log(-). 

Both inequalities play a key role in the next two chapters, which address inequalities on the 
two-dimensional sphere § 2 . It is worth noting that Ghigi's inequality relies on the Prekopa- 
Leindler principle, which itself is another manifestation of a mass transport context. One can 
therefore infer that the approach of Part IV can and should be made more readily applicable 
to critical Moser-type inequalities. 

In Chapter 17, we establish the Moser-Trudinger inequality, which states that for a > 1, 
the functional 

(0.35) J a (u):=a \Vu\ 2 duj + 2 udw -In / e 2u duj 

Js 2 Js 2 Js 2 

is bounded below on the Sobolev space i? 1,2 (§ 2 ), where here du := j^sin9d9 A dip denotes 
Lebesgue measure on the unit sphere, normalized so that J s2 du> = 1. We also give a proof of 
Onofri's inequality which states that the infimum of J a on H 1,2 (S 2 ) is actually equal to zero 
for all a > 1, and that 

(0.36) inf{J x (it); u G tf 1 * 2 ^ 2 )} = inf{J x (tt); ueM} = 0, 

where M is the submanifold M = {u E if 1,2 (§ 2 ); J g2 e 2u xduj = 0}. Note that this inequality, 
once applied to axially symmetric functions, leads to the following counterpart of (|0.33[) 

(0.37) inf I a (g) = inf IJg) = if a > 1. 

gtH 1 (-1,1) 

In Chapter 18, we address results of T. Aubin asserting that once restricted to the submanifold 
Ai, the functional J a then remains bounded below (and coercive) for smaller values of a, which 
was later conjectured by A. Chang and P. Yang to be equal to |. We conclude the latest 
developments on this conjecture, including a proof that 

(0.38) mi{J a (u); u E M} = ifa>§. 

The conjecture remains open for 1/2 < a < 2/3. 
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We have tried to make this monograph as self-contained as possible. That was not possible 
though, when dealing with the applications such as in Chapters 4, 7 and 10. We do however 
give enough references for the missing proofs. 

The rapid development of this area and the variety of applications forced us to be quite 
selective. We mostly concentrate on certain recent advances not covered in the classical books 
such as the one by R. A. Adams [3] and V. G. Maz'ya |167j . Our choices reflect our taste 
and what we know -of course- but also our perceptions of what are the most fundamental 
functional inequalities, the novel methods and ideas, those that are minimally ad-hoc, as well 
as the ones we found useful in our work. It is however evident that this compendium is far from 
being an exhaustive account of this continuously and rapidly evolving line of research. One 
example that comes to mind are inequalities obtained by interpolating between the Hardy and 
the Trudinger-Moser inequalities. One then gets the singular Moser-type inequalities, which 
states that for some C = C (n, |ft|) > 0, one has for any u S W ' n (ft) with J Q \Vu\ n dx < 1, 

r exp(/?M^) 
(0.39) / Ux<C , 

1 n 

for any a € [0, n) , < (3 < (l — ^) noj^z\ , where w„_i = is the area of the surface of 

the unit n-dimensional ball. See for instance [H [9]). 

The recent developments on these inequalities could have easily constituted a Part VII 
for this book, but we had to stop somewhere and this is where we stop. 
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